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\Q ABSTRACT 

The results of analysis of HST /—band imaging of a sample of 41 radio galaxies span- 
ning three orders of magnitude in radio luminosity at redshift z ~ 0.5 are presented. 
These results represent the first stage in a coordinated programme to study the con- 
nection between radio luminosity and host-galaxy properties, black-hole masses and 
cluster environments in radio galaxies spanning a wide range in radio luminosity over 
a restricted range in redshift. The full sample is comprised of objects drawn from four 
complete, low-frequency selected radio samples with progressively fainter flux-density 
Qh| limits (3CRR, 6CE, 7CRS and the new TexOx-1000 sample). Modelling of the HST 

q ■ imaging data shows that the host galaxies have surface-brightness distributions con- 

5— i ' sistent with those expected for classic ellipticals (Sersic parameter, ~ 0.25), with j3 

\ in the range 0.17 < (3 < 0.30, and a mean of < (3 >= 0.23 ± 0.01. The luminosities 

& ■ of the host galaxies are found to be comparable with those of galaxies drawn from 

the bright end of the local cluster galaxy luminosity function, spanning the range 
0.7L* < L < 10L*, with a mean of 3.2 ± 0.3L*, after correcting for the effects of pas- 
sive evolution. In addition, the radio galaxies are shown to follow a Kormendy ([i e — r e ) 
relation indistinguishable from that of both powerful low-redshift radio galaxies and 
local Abell brightest cluster galaxies. Combining our new results with those in the 
literature it is found that the scalelengths and Kormendy relations of 3C-class radio 
galaxies do not vary significantly over the redshift range 0.0 < z < 0.8, providing 
no evidence for dynamical evolution of this class of host galaxy within this redshift 
interval. Converting the host-galaxy luminosities into black-hole mass estimates, using 
the local Mf,h — Mt u i ge correlation, predicts that the radio galaxies harbour central 
black holes with masses in the range 10 81 Mq < M^h < 1O 9,5 M , with a geometric 
mean of < Mf,h >= 10 8 87±0 04 M Q . Finally, a significant (~ 3a) correlation is found 
between black-hole mass and 151-MHz radio luminosity for those objects in the sam- 
ple with either high-excitation nuclear spectra (HEG) or classical double (CD) radio 
structures. 
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1 INTRODUCTION 



The underlying mechanisms which govern the large range 
in low-frequency radio luminosity observed in radio-loud ac- 
* Email: rjm@roe.ac.uk tive galactic nuclei (AGN) are currently unknown. Over the 
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previous two decades a wealth of observational effort has 
been invested in studying the redshift regime < z < 1, 
investigating whether the dominant factors are linked to 
the Mpc-scale cluster environments (eg. Prestage & Pea- 
cock 1988; Hill & Lilly 1991; Ellingson, Yee & Green 1991; 
Wold et al. 2000; Best 2000) or the kpc-scale properties of 
the host galaxies (e.g. Smith & Heckman 1989; Best, Lon- 
gair & Rottgering 1998; McLure & Dunlop 2000; Dunlop et 
al. 2003; Zirm, Dickinson & Dey 2003). 

One of the incontrovertible observational facts which 
has emerged from these studies, at least at z < 1, is that 
the host galaxies of all powerful radio-loud AGN are mas- 
sive L > L* ellipticals (eg. Taylor et al. 1996; McLure et al. 
1999; Dunlop et al. 2003). The apparent uniformity of radio- 
loud AGN host galaxies has taken on added importance 
over the last few years, following the discovery in nearby 
(distance < 150 Mpc) inactive galaxies that a reasonably 
accurate estimate (AMjj, ~ 0.3 dex) of the central black- 
hole mass can be obtained via its correlation with the mass 
of the host spheroidal component (Kormendy & Richstone 
1995; Magorrian et al. 1998; Gebhardt et al. 2000; Ferrarese 
& Merritt 2000; McLure & Dunlop 2002; Marconi & Hunt 
2003; Tremaine et al. 2002). Moreover, recent progress has 
also indicated that a similarly accurate black-hole mass esti- 
mate (AMbh — 0.4 dex) can be obtained for broad-line AGN 
using emission-line widths to derive the virial mass estimate 
(eg. Kaspi et al. 2000; McLure & Dunlop 2002; McLure & 
Jarvis 2002; Vestergaard 2002). Consequently, a large body 
of work has appeared in the recent literature investigating 
the possible link between radio luminosity and black-hole 
mass in radio-loud AGN (eg. Laor 2000; Lacy et al. 2001; 
McLure & Dunlop 2001a; McLure & Dunlop 2002; Bettoni 
et al. 2003; Dunlop et al. 2003). 

Unfortunately, observational studies have traditionally 
been subject to a degeneracy between radio luminosity and 
redshift produced as a by-product of flux-limited radio sam- 
ples. In order to study the properties of radio-loud AGN 
separated by a large dynamic range in radio luminosity, it 
has previously been necessary to select samples consisting 
of objects covering a wide range of redshifts. This has led 
to difficulties in interpreting the data due to the complica- 
tion of potentially significant evolutionary effects. However, 
by selecting our sample of objects from four complete, low- 
frequency selected radio samples with successively fainter 
flux-density limits, it has been possible to construct a sample 
of radio galaxies which spans three decades in radio luminos- 
ity at a virtually constant cosmic epoch (0.4 < z < 0.6). In 
this respect the motivation for this study is identical to that 
pursued by Hill & Lilly (1991), who investigated the clus- 
ter environments of a sample of 45 radio galaxies spanning a 
similar range in radio luminosities at z ~ 0.5. By successfully 
disentangling the effects of luminosity and redshift for the 
first time, Hill & Lilly (1991) were able to demonstrate that 
z ~ 0.5 marks an apparent epoch-dependent change in clus- 
ter environment, with radio galaxies tending to inhabit sig- 
nificantly richer cluster environments than their low-redshift 
counterparts. Furthermore, the wide range of radio luminosi- 
ties in the Hill & Lilly sample enabled them to determine 
that extended radio luminosity did not appear to be a strong 
function of cluster richness. 

Developments over the intervening decade have pro- 
vided our new study with two crucial advantages over the 



previous work of Hill & Lilly (1991). Firstly, the availabil- 
ity of the 7CRS (Willott et al. 2003; Lacy et al. 1999) 
and TexOx-1000 (Hill & Rawlings 2003) radio samples has 
now allowed us to compile a sample of 41 radio galaxies at 
z ~ 0.5, which is complete over the full three-decades range 
in radio power. Secondly, with the high-resolution imaging 
provided by the Hubble Space Telescope (HST), and the im- 
provement in follow-up observations available with ground- 
based telescopes, we are now in a position to obtain a wide 
range of high quality data-sets for our new sample which 
were not obtainable only ten years ago. Consequently, our 
new sample of z ~ 0.5 radio galaxies is the basis for an am- 
bitious project, which will for the first time systematically 
investigate the connection between the radio luminosity of 
radio galaxies, the properties of their host galaxies, their lo- 
cation on the fundamental plane, the richness and mass of 
their cluster environments and the masses of their central 
black holes. 

The results presented in this paper represent the first 
stage of this wide-ranging project, and concentrate solely 
on the analysis of our deep /—band HST WFPC2 imaging 
data. The results of our extensive follow-up observations of 
the sample are deferred to a series of future papers. 

The structure of the paper is as follows, in Section 2 
the details of the radio-galaxy sample are described, while 
in Section 3 the HST observations and data reduction are 
discussed. In Section 4 the modelling of the HST imaging 
is described, with the main results presented and compared 
to literature results on low-redshift radio galaxies in Sec- 
tion 5. In Section 6 the host-galaxy properties of our 3CRR 
sub-sample are compared to literature results to explore the 
evidence for dynamical evolution of the most powerful radio 
galaxies within the redshift range 0.0 < z < 0.8. In Sec- 
tion 7 the properties of the radio-galaxy hosts are compared 
to those of local brightest cluster galaxies. In Section 8 the 
masses of the radio galaxies' central black holes are esti- 
mated, and in Section 9 the relationship between black-hole 
mass and radio luminosity is investigated. The main conclu- 
sions are summarized in Section 10. Unless otherwise stated, 
throughout this paper the following cosmology is assumed: 
H = 70 km s _1 Mpc _1 , fl m = 0.3, Q A = 0.7. 



2 THE SAMPLE 

The full 41-object z ~ 0.5 radio-galaxy sample (hereafter 
the ZP5 sample) consists of all the narrow-line radio galax- 
ies in the redshift interval 0.4 < z < 0.6 from four, com- 
plete, low-frequency selected radio surveys; 3CRR (Laing, 
Riley & Longair 1983), 6CE (Eales et al. 1997; Rawlings, 
Eales & Lacy 2001), 7CRS (Willott et al. 2003; Lacy et al. 
1999) and TexOx-1000; hereafter TOOT (Hill & Rawlings 
2003). The locations of all four radio surveys on the radio 
luminosity — redshift plane are shown in Fig 1. Details of 
the ZP5 sample can be found in Table 1. Additionally, Table 
2 lists those objects which satisfied the redshift and radio lu- 
minosity criteria for the ZP5 sample, but were excluded due 
to their broad-emission line spectra. 

The TOOT objects were selected from a preliminary 
version of the survey. They are therefore a subset of all the 
TOOT sources meeting the radio flux density and sky area 
selection criteria. Broad-line objects were excluded, in the 
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Table 1. The full ZP5 sample. Column one lists the radio-galaxy names, while columns two and three list the J2000 source coordinates. 
Column 4 lists the radio-galaxy rcdshifts and column five lists the logarithm of the 151-MHz luminosities in units of WHz _1 sr _1 . Column 
6 classifies the radio structures as either Classical Double (CD), Fat Double (FD), Jet (J) or Compact (CM), the same classifications 
as used by Owen & Laing (1989). The radio structure classifications are based on consistent 1.4 GHz VLA A-array data (1.5 arcsec 
resolution) for the 6CE, 7CRS and TOOT sub-samples. For the 3CRR sub-sample, for which structural classifications are well established, 
the classifications were based on the high-quality maps available in the literature (at a variety of frequencies/resolutions). Column 7 
describes the optical nuclear spectra in terms of the presence of high-excitation narrow emission lines (HEG) or not (LEG), using the 
classification scheme of Jackson & Rawlings (1997). Those objects for which the available optical spectra are of insufficient quality to 
allow a reliable HEG/LEG determination are listed with a '?'. To avoid confusion we include here in parenthesis the old names of the 
7CRS objects which have been adopted in previous papers: 7C0213+3418 (5C6.63), 7C0219+3419 (5C6.201), 7C0219+3423 (5C6.214), 
7C0220+2952 (5C6.233), 7C0223+3415 (5C6.279), 7C0810+2650 (5C7.7). 
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case of TOOT 1252+3310 after HST imaging and during fur- 
ther spectroscopic follow-up. There is some possibility that 
this subset is biased in the sense that the optically faintest 
sources may be have been preferentially excluded, although 
the level of bias will only be clear once analysis of the full 
survey is complete. 

The original decision to select the radio-galaxy sample 



in the redshift interval 0.4 < z < 0.6 was motivated by 
several factors. Principal among these was that at z ~ 0.5 
it becomes possible to get a full three decades of dynamic 
range in radio luminosity, including objects which are close 
to the most radio luminous objects which exist. Secondly, at 
z ~ 0.5 it is possible to construct sub-samples, divided by ra- 
dio luminosity, which contain sufficient numbers of objects to 
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Figure 1. 151-MHz radio luminosity versus redshift for the narrow-line radio galaxy component of the 3CRR, 6CE and 7CRS samples. 
The members of the ZP5 sample are highlighted as enlarged symbols in the redshift interval 0.4 < z < 0.6, including those drawn from 
the new TexOx-1000 (TOOT) survey. 



Table 2. List of sources which met the redshift and radio-luminosity criteria for inclusion in the ZP5 sample, but were excluded due to 
the presence of broad emission lines. In column 2 the sources are classified as cither quasar (Q) or weak quasar (WQ) according to the 
classification scheme of Jackson & Rawlings (1997). We note here that the excluded quasars are confined to the two most radio luminous 
sub-samples, in agreement with the findings of Willott et al. (2000) that the quasar fraction in complete, low-frequency selected radio 
samples is an increasing function of radio luminosity. The implications of the possible quasar cores in the ZP5 sample (see Table 3) will 
be fully discussed in a forthcoming paper (Mitchell ct al., in prep). 



Source 
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Comment 


3C 47 


Q 




3C 147 


Q 




3C 215 


Q 




3C 275.1 


Q 




3C 334 


Q 




3C 455 


WQ 




6C 1220+3723 


Q 


B1950.0 name - as in Rawlings ct al. (2001) 


5C 7.118 


WQ 


B1950.0 name - as in Willott et al. (2003) 


TOOT 1252+3310 


WQ 


Broad H/3 detected in follow-up spectroscopy (Willott ct al., in prep) 



reliably determine the mean values and variances of relevant 
parameters, allowing any correlations with radio luminosity 
to be investigated. Thirdly, the chosen redshift range makes 
it possible to obtain the deep, high-resolution images nec- 
essary to derive accurate host-galaxy parameters, without 
requiring prohibitive amounts of HST time, while the upper 
limit of z = 0.6 ensures that with the WFPC2 filters we are 
still able to sample light longward of the 4000A break. Fi- 
nally, redshift z ~ 0.6 represents the maximum redshift for 
which it is practicable on 4m-class telescopes to obtain accu- 
rate measurements of the host-galaxy stellar-velocity disper- 
sions, necessary to investigate the radio-galaxy fundamental 
plane properties (Willott et al., in prep). 



3 OBSERVATIONS AND DATA REDUCTION 

The observations for 39/41 of the objects in this project 
were carried out in the Cycle 10 program GO#9045 with 
the Hubble Space Telescope. We used the WFPC2 detec- 
tor with the F785LP filter. When combined with the sys- 
tem efficiency this filter has a peak wavelength of 8700 A 
and high throughput from 7900 A- 9300 A. This filter was 
selected because we wish to be sampling light above the 
4000 A break, where there is less likely to be contamination 
of the old-stellar population by more recently formed stars 
or AGN-related emission. Further restrictions on the filter 
choice came from our desire to be free from emission lines 
as much as possible. The [O ill] A5007 line only enters the 
filter bandpass at redshifts z > 0.58 and Ha leaves the filter 
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bandpass at redshifts z > 0.42. Since our sample spans the 
redshift range 0.4 < z < 0.6 this ensures very few sources 
will have significant emission-line contamination with this 
filter. 

The radio galaxies were placed in the WF3 chip so that 
information on their arcminute-scale environments could be 
probed with the other chips on WFPC2. The orientation 
of the detector on the sky was arranged so that diffraction 
spikes and blooming from bright stars did not affect the 
array at the locations of the radio galaxies. The observations 
of each radio galaxy were made during a single orbit and 
the total exposure time of 2000 seconds was split into 4 
integrations of 500s each to facilitate the removal of cosmic 
rays. The target was dithered 5.5 WF3 pixels in x and y in 
between integrations to improve the sampling of the PSF. 

Calibrated images were retrieved from the HST archive. 
Registering and combining the 4 images per target were 
performed using the drizzle routines in the iraf package 
(Fruchter & Hook 2002). The final reduced images for all 
four CCDs were sampled on a grid of pixels with size 0.0498 
arcsec. Unsaturated stars typically have FWHM values of 3 
pixels corresponding to 0.15 arcsec. 

Existing observations were used in the analysis of the 
two remaining objects in the sample; 3C16 and 3C295. In the 
case of 3C16, an existing HST WFPC2 observation using the 
PC chip and the F702W filter was obtained from the HST 
archive (GO#6675). The pixel size of the PC chip (0.045 
arcsec/pixel) is well matched to that of the final drizzled 
images of the rest of the sample, and because 3C16 is at the 
low-redshift end of the sample (z = 0.405), the F702W filter 
is still sufficiently red to sample light longward of the 4000A 
break. Although an archival HST observation of 3C295 was 
also available, because of the placement of the radio galaxy 
close to the edge of the PC chip, it was deemed unsuitable 
for the host-galaxy analysis described in Section 4. Fortu- 
nately we were provided with a deep J— band image of the 
3C295 cluster obtained on the 2.5m Nordic Optical Tele- 
scope for the weak lensing analysis of Wold et al. (2002). 
The final image comprises 10,000 seconds of integration and 
has a FWHM of 0.8 arcsecs. Although this is clearly lower 
spatial resolution than the HST data, this does not have any 
influence on the accurate determination of the host-galaxy 
properties because 3C295 is the largest, most luminous tar- 
get in the whole sample. 



4 HOST-GALAXY ANALYSIS 

Given the high resolution and depth of the HST data avail- 
able for this study it was decided that the determination 
of the host-galaxy parameters should be carried out via 
full two-dimensional modelling of the images. The two- 
dimensional modelling technique adopted is identical to that 
which has previously been applied in the analysis of both ra- 
dio galaxies and quasar host galaxies (eg. McLure et al. 1999; 
McLure & Dunlop 2000; Dunlop et al. 2003). Full details of 
this modelling technique have been previously published in 
McLure, Dunlop & Kukula (2000) and so only a brief outline 
will be provided here. 

The first stage in the process is the production of a 
two-dimensional model host galaxy which has a surface- 



brightness distribution governed by an axially-symmetric 
Sersic (1968) profile of the form: 



E(r) = So exp 



(- 



(1) 



where So is the central surface brightness, ro is the charac- 
teristic scalelength and the ft parameter controls the shape 
of the radial profile. Under the adopted form of the Sersic 
profile a value of /? = 0.25 corresponds to a classical de 
Vaucouleurs r 1 ^ 4 law (de Vaucouleurs 1953) and a value 
of f3 = 1.0 corresponds to an exponential Freeman (1970) 
disc model. This model galaxy is subsequently convolved 
with the HST point spread function in order to produce 
a synthetic host-galaxy image. A chi-square comparison of 
the data and model is then performed on a pixel-by-pixel 
basis, in order to determine the goodness of fit of the ini- 
tial model. A chi-square minimization process is then per- 
formed to determine the best-fitting value of the six param- 
eters which describe the model (Eo,ro,/3, axial ratio, posi- 
tion angle and nuclear point-source luminosity). Although 
the model galaxy axial ratio and position angle are included 
as free parameters during the fitting process, they are held 
constant as a function of radius. During the model fitting 
process it was not assumed a priori that the radio galaxies 
would not possess an unresolved nuclear point source. 



4.1 Point spread function 

For the analysis of the HST imaging data we utilized the 
model WFPC2 point spread functions (PSF) produced by 
the TinyTim software package (Tiny Tim v6.0, Krist 2001). 
The great advantage of using the TinyTim model PSFs over 
empirical WFPC2 PSFs is that they have infinite signal- 
to-noise, and can be produced to match the (x, y) coordi- 
nates of each radio galaxy on the WF3 chip. Although the 
model PSFs produced by TinyTim do not reproduce the 
large-scale halo of scattered light of the empirical WFPC2 
PSF at radii of r > 1.5", this was not considered to be a 
serious problem for the current analysis because the lumi- 
nosity of the any nuclear point sources were anticipated to 
be only a small fraction of the host-galaxy luminosity. This 
assumption proved to be justified, with the vast majority 
of the radio galaxies proving to possess a negligible nuclear 
point-source component (the median nuclearihost flux ratio 
for the ZP5 sample is 0.03). 

In order to match the final reduced radio-galaxy im- 
ages as accurately as possible, model PSFs were produced 
at the two dithered positions of each radio-galaxy observa- 
tion separately. The two model PSFs were then processed 
by the drizzle routine in an identical fashion to the radio- 
galaxy images to produce the final model PSF used during 
the modelling process. 



5 RESULTS 

In this section the results of the two-dimensional modelling 
of the host galaxies of the ZP5 sample are presented and 
discussed. A list of the principal host-galaxy parameters de- 
rived from the model fits for each individual object can be 
found in Table 3, while the average values for parameters 
of particular interest for the four radio-galaxy sub-samples 
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Figure 2. Examples of the two-dimensional modelling analysis of the HST data for one object from each of the four sub-samples. The 
left-hand panels show cut-outs of the reduced HST WFPC2 images centred on the radio galaxies. The middle panels show the best-fitting 
host-galaxy models. The right-hand panels show the residual images produced by subtracting the best-fitting models from the reduced 
radio-galaxy images. All panels are 12.75" X 12.75". 
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and the full ZP5 sample can be found in Table 4. A graph- 
ical representation of the modelling process is shown in Fig 
2 which, for one radio galaxy from each sub-sample, shows 
postage-stamp images of the reduced HST data, the best- 
fitting two-dimensional models and the residual images after 
subtraction of the best-fitting models. 



5.1 Morphology 

The mean value of the Sersic /3 parameter for the full ZP5 
sample is < (3 >= 0.23 ±0.01, consistent with that expected 
for classical elliptical galaxies (the typical uncertainty in (3 
is ±0.05). The average (3 values displayed in Table 4 for 
the four radio-galaxy sub-samples show no evidence for any 
trend between the (3 parameter and radio luminosity. 

In Fig 3 we show a histogram of the derived (3 param- 
eters for the full ZP5 sample. Although Fig 3 confirms that 
the morphologies of the ZP5 sample radio galaxies are con- 
sistent with the classical j3 = 0.25 value, we note here that 
there are a substantial number of objects with [3 < 0.25. 
The possible implications of this in light of the recent discov- 
ery of a correlation between (3 and central black-hole mass 
(Graham et al. 2001; Erwin et al. 2003) will be discussed in 
Section 8.1. 

The distribution of (3 parameters derived for the ZP5 
sample is in good agreement with that determined by Dun- 
lop et al. (2003) for their sample of twenty 3C-class radio 
galaxies and radio-loud quasars at z ~ 0.2. The mean (3 
parameter determined by Dunlop et al. for their radio-loud 
AGN sample was < j3 >= 0.24 ± 0.01. Unfortunately there 
are no comparable samples in the literature of radio galaxies 
at higher redshift, for which Sersic f3 parameters have been 
derived, with which to compare the morphological results 
derived for the present sample. 



5.1.1 Model- subtracted images 

In the final column of Table 3 brief comments are listed 
describing the structure revealed by subtracting the best- 
fitting axially-symmetric host-galaxy model from the radio- 
galaxy images (see Fig 2 for examples). Objects for which 
the best-fitting galaxy model produces a clean subtrac- 
tion, with little or no residual structure present in the sub- 
tracted images, are described as being undisturbed (UD). 
These objects comprise one third of the full ZP5 sam- 
ple. Objects which show clear signs of interaction are de- 
scribed as being mergers (M), with those cases where 
the secondary nucleus/companion object has a luminosity 
> 10% of that of the radio-galaxy being described as ma- 
jor mergers (MM). Objects for which the modelling indi- 
cates the presence of an unresolved point-source nucleus 
with luminosity > 10% of that of the radio galaxy (eg. 
TOOT 1303±3334, 7C 1731±6638) are described as possi- 
ble quasar cores (QC?). Finally, those objects whose model- 
subtracted images reveal symmetric, central positive flux 
residuals are described as having possible disc or aligned 
emission (DA?). 
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Figure 3. Histogram of the derived Sersic parameters for the 
full ZP5 sample. The dotted vertical lines indicate the locations 
of f3 = 0.25 and (3 = 1.0, the expected values for classic de Vau- 
coulcurs r 1 / 4 and Freeman exponential disk profiles respectively 
(de Vaucouleurs 1953; Freeman 1970). 



The information provided by the model-subtracted im- 
ages is not pursued further within this paper, which con- 
centrates on the properties of the dominant smooth stellar 
component of the radio-galaxy hosts. However, this infor- 
mation will be fully exploited by a forthcoming paper which 
explores the link between radio luminosity /structure and the 
kpc-scale environments of the radio galaxies, the incidence 
and magnitude of interaction/merger activity, aligned emis- 
sion components and the properties of those radio galaxies 
with possible quasar cores (Mitchell et al., in prep). 



5.2 Host-galaxy luminosity 

The absolute magnitudes of the radio galaxies have been cal- 
culated under the assumption that the spectral energy distri- 
butions of the hosts can be well represented by the passively 
evolving stellar populations typical of massive early-type 
galaxies. This decision is based on several lines of evidence, 
including the morphological results presented in the previ- 
ous section, the R — K colour versus redshift distribution 
of 7CRS objects (Willott, Rawlings & Blundell 2001) and 
the off-nuclear spectroscopy of z ~ 0.2 quasar host galaxies 
(Nolan et al. 2000). Furthermore, the initial results from our 
follow-up spectroscopic observations, designed to determine 
the host-galaxy stellar- velocity dispersions, indicate that the 
host-galaxy spectra are consistent with mature stellar pop- 
ulations (Willott et al., in prep). 

The absolute magnitudes for each host galaxy pre- 
sented in Table 3 are based on the integrated instrumen- 
tal (F785LP) magnitudes returned by the modelling code. 
These instrumental magnitudes are then converted to appar- 
ent 7— Cousins magnitudes using the filter transformations 
and colour equations detailed in Holtzman et al. (1995). The 
apparent colours required for this transformation were cal- 
culated from the 1995 version of Bruzual & Chariot single- 
burst elliptical galaxy stellar-population models (Bruzual & 
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Table 3. Host-galaxy parameters derived from the two-dimensional modelling of the HST WFPC2 data. Column one lists the radio- 
galaxy names. Column 2 lists the effective radii in kpc. Column 3 lists the mean R— band surface brightness within the effective radii 
in mag.arcsec -2 . The surface brightnesses have been corrected for (1 + z) 4 cosmological dimming, and then K-corrected and corrected 
for passive evolution as described in the text. Column four lists the axial ratios and column five lists the real optical position angles of 
the radio galaxies in degrees, transformed from the fitted on-chip position angles. Column six lists the Sersic (3 parameters. Column 7 
lists the integrated absolute R— band host luminosities, including K-correction and correction for passive evolution as described in the 
text. Column 8 lists the host-galaxy luminosities in units of L* , where L* is taken to be Mr = —22.07 (Yagi et al. 2002). Column 9 
lists the ratio of apparent nucleanhost luminosity. Column 10 lists the estimated black-hole mass for each object based on its fitted 
bulge luminosity as described in Section 8. Column 11 highlights properties of the radio galaxies revealed by the model-subtracted 
images: UD=undisturbed, M=merger, MM=major merger, DA?=possible disc/aligned emission, QC?=possible quasar core (see text for 
discussion). 
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Chariot 1993), under the assumption of a z = 5.0 forma- 
tion redshift. Following correction for cosmological dimming 
within our chosen cosmology, evolution and K— corrections 
were then applied using the same Bruzual & Chariot mod- 
els. Galactic extinction corrections were then applied to each 
object (Schlegel et al. 1998), before a final conversion to ab- 
solute R— Cousins magnitudes was made, assuming a z — 
colour of R — I = 0.7 (Bruzual & Chariot 1993), in order 



to facilitate a straightforward comparison with literature 
results on radio-galaxy properties. The surface-brightness 
results presented in Tables 3 & 4 have been corrected for 
(1 + z) 4 cosmological dimming, and then K-corrected and 
corrected for passive evolution in the same fashion as the ab- 
solute magnitudes. The host-galaxy luminosities presented 
in Tables 3 & 4 and Fig 4 demonstrate that the host galax- 
ies of all four sub-samples are luminous, with a mean inte- 
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Table 4. Average values of the host-galaxy parameters derived from the two-dimensional modelling for the four radio-galaxy sub-samples 
and the ZP5 sample as a whole. For each parameter the mean value for the relevant sample is listed with its associated standard error, 
and the median value is listed immediately below. Column 2 lists the Sersic /3 parameters. Column 3 lists the axial ratios, while column 
4 lists the effective scalelengths in kpc. Column 5 lists the mean R— band surface brightness within the effective radius in mag.arcsec -2 . 
Column 6 lists the absolute integrated R— band magnitudes, while column 7 lists the host-galaxy luminosities in units of L* . 
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Figure 4. Histogram of integrated absolute R— band host mag- 
nitudes for the full ZP5 sample. The vertical dotted line indicates 
the location of L* for local cluster galaxies: Mr = —22.07 (Yagi 
ct al. 2002). 

grated R— band absolute magnitude for the full ZP5 sample 
of < Mr >— —23.20 ± 0.09. In fact, these luminosities are 
directly comparable to the luminosities of galaxies drawn 
from the bright end of the local cluster-galaxy luminosity 
function. A comparison with the recent study of Yagi et al. 
(2002) shows that the mean host-galaxy luminosity of the 
full ZP5 sample is 3.2 ± 0.3L*, where the Yagi et al. deter- 
mination of L* for local cluster galaxies is Mr = —22.07 in 
our adopted cosmology. 

The mean host-galaxy luminosities presented in Table 
4 also allow a comparison to be made with the host-galaxy 
luminosity differences between the 3C, 6C and 7C samples 
determined from the K — z diagram (Willott et al. 2003; 
Eales et al. 1997). In the if -band, Willott et al. find a 
mean luminosity difference of 0.33 magnitudes between the 
3C and 6C samples, and a corresponding difference of 0.55 
magnitudes between the 3C and 7C samples. The equivalent 



figures from the analysis of our HST imaging are 0.34±0.16 
magnitudes and 0.76 ± 0.25 magnitudes in the R— band, in 
good agreement with the K— band results. 

The rough scaling between extended radio luminosity 
and host-galaxy luminosity apparent from both the K — z 
relation and the new results presented here, suggests a cor- 
relation exists between radio luminosity and black-hole mass 
within the 3CRR, 6CE and 7CRS samples. This issue is dis- 
cussed further in Section 9. 



5.3 Host-galaxy scalelengths 

The best-fitting scalelength from the two-dimensional mod- 
elling of each individual member of the ZP5 sample is listed 
in Table 3, while the mean and median values for the four 
radio-galaxy sub-samples and the full ZP5 sample are listed 
in Table 4. 

The full range of scalelengths displayed by the objects 
in the ZP5 sample is 5 < r e < 31 kpc, with a mean value 
of 13.9 ± 0.9 kpc. Scalelengths of this size clearly place the 
radio galaxies in the regime of brightest cluster galaxies (eg. 
Graham et al. 1996; Hoessel, Oergerle & Schneider 1987) 
and a comparison of the ZP5 host galaxies with the proper- 
ties of brightest cluster galaxies (BCGs) is pursued further 
in Section 7. Furthermore, in Section 6 the scalelength re- 
sults for the 3CRR sub-sample are combined with existing 
literature results on 3C-class galaxies to search for any ev- 
idence for dynamical evolution amongst the most powerful 
radio galaxies within the redshift interval 0.0 < z < 0.8. 

In common with the results regarding the ZP5 host- 
galaxy luminosities, an inspection of Table 4 reveals an ap- 
parent correlation between mean scalelength and radio lumi- 
nosity within the 3CRR, 6CE and 7CRS sub-samples. This 
is entirely as expected given the correlation between host- 
galaxy and radio luminosity discussed above, and the well 
known correlation between host luminosity and scalelength 
among early-type galaxies (e.g. r e oc £°'o S 3 t ±0 ' 03 ; Bernardi et 
al. 2003) 
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Figure 5. The top panel shows the best-fitting Kormendy rela- 
tion for the full ZP5 sample, before the application of any evolu- 
tionary corrections. The middle panel shows the same ZP5 data 
with the addition of the 73 low-rcdshift radio galaxies from Bet- 
toni et al. (2001). The best-fitting fixed-slope Kormendy relations 
to the ZP5 and Bettoni et al. radio galaxies are shown as the 
solid and dashed lines respectively (see text for discussion). The 
bottom panel shows the best-fitting Kormendy relation to the 
combined sample of the ZP5 and Bettoni et al. radio galaxies. In 
the bottom panel the ZP5 host galaxies have been corrected for 
passive evolution of their stellar populations as described in the 
text. 



5.4 The Kormendy relation 

Early-type galaxies are known to exist on a two-dimensional 
manifold (fundamental plane) in the three-dimensional pa- 
rameter space defined by effective scalelength, the mean 
surface brightness within the effective scalelength and the 
central stellar-velocity dispersion (eg. Dressier, Lynden-Bell 
& Burstein 1987; Djorgovski & Davis 1987). At present we 
are engaged in obtaining stellar- velocity dispersion measure- 
ments for the ZP5 sample to investigate their location on the 
fundamental plane, and to estimate the masses of their cen- 
tral black holes via the Mth — o correlation (Willott et al., in 
prep) . However, armed with just the scalelength and surface 
brightness parameters for each host galaxy, it is possible to 
examine the photometric projection of the full fundamental 
plane, the so-called Kormendy or fi e — r e relation (Kormendy 
1977). 

The Kormendy relation for the ZP5 sample is shown in 
the top panel of Fig 5. In this panel the < fi > e values of the 
ZP5 objects have been K-corrected and corrected for (1 + z) 4 
surface-brightness dimming, but have not been corrected for 
passive evolution. Using the iterative \ 2 FITEXY routine, 
which takes account of errors in both parameters (Press et 
al. 1992) , the best-fitting form of the Kormendy relation for 
the non-evolutionary corrected ZP5 sample is: 



< n > e = 3.72(±0.51) logr e + 16.08(±0.59) 



(2) 



In the middle panel of Fig 5 the same ZP5 data are repeated, 
along with the 73 low-redshift radio galaxies from the funda- 
mental plane study of Bettoni et al. (2001). The best-fitting 
relation to Bettoni et al. sample is: 



< ii > e = 3.32(±0.14) logr e + 16.97(±0.15) 



(3) 



Given that the Bettoni et al. sample is virtually at redshift 
zero, the vertical offset between the Kormendy relation of 
the Bettoni et al. sample and that of the ZP5 sample is a 
direct indication of the amount of luminosity evolution be- 
tween z — and z — 0.5. If the two samples are fitted with 
a Kormendy relation of fixed slope 3.52 (intermediate to 
their two independent fits) then the best-fitting intercepts 
are 16.30(±0.07) and 16.77(±0.05), for the ZP5 and Bettoni 
et al. samples respectively. This directly implies 0.47 ± 0.09 
magnitudes of R— band luminosity evolution between z=0.5 
and z=0. This is in excellent agreement with the evolution- 
ary corrections derived from the Zf or = 5 Bruzual & Chariot 
model, which predicts a mean value of 0.42 ± 0.01 magni- 
tudes of passive luminosity evolution for the ZP5 sample. 

In the bottom panel of Fig 5 the low-redshift radio 
galaxies of Bettoni et al. are again plotted, together with 
the now evolutionary corrected ZP5 objects. The solid line 
shows the best-fitting relation to the combined sample which 
has the form: 



< fi > e = 3.35(±0.13) logr e + 16.93(±0.14) 



(4) 



which is in excellent agreement with the recent determina- 
tion of the Kormendy relation slope (3.33 ± 0.09) for 9000 
z ~ early-type galaxies drawn from the Sloan Digital Sky 
Survey (SDSS) by Bernardi et al. (2003). 
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Figure 6. Panel A shows the Kormendy relation for the four 3C-class radio-galaxy samples described in the text. The surface brightnesses 
have been K-corrected and corrected for passive evolution as described in Section 5.2. The solid line is the best-fitting relation to the 
combined sample and has a slope of 3.78 ± 0.30. Panel B shows the evolution of the mean absolute R— band magnitudes (with standard 
error bars) of the four 3C-class samples with redshift, where the absolute magnitudes have been K-corrected, but have not been corrected 
for passive evolution. The solid line shows the expected evolution of the mean absolute magnitude with redshift due to the passive 
evolution of their stellar populations, using the Bruzual & Chariot single burst elliptical models with Zf or = 5 and assuming a redshift 
zero absolute magnitude of Mr = —23.35. 



6 EVOLUTION OF THE MOST POWERFUL 
RADIO GALAXIES 

In this section we focus on our new results for the 3CRR 
ZP5 sub-sample, combining them with previous work in the 
literature, to re-examine the issue of the evolution of 3C- 
class radio galaxies from z ~ 0.8 through to the present day. 
It has been known for some twenty years now that powerful 
radio galaxies follow a tight correlation between their ap- 
parent K— band magnitudes and redshift (Lilly & Longair 
1984). The simplest interpretation of this K — z relation con- 
sistent with the data assumes that the radio galaxies were 
formed at high redshift (zf or > 3), either in a monolithic 
collapse or via hierarchical merging, and have been evolving 
passively thereafter. 

However, this interpretation of the K — z relation was 
challenged by Best, Longair & Rottgering (1998) who argued 
that the apparent passive appearance of the K — z relation 
was in fact a conspiracy. Based on evidence in the literature 
that points to higher redshift 3C radio galaxies (z > 0.5) 
inhabiting significantly richer cluster environments (Hill & 
Lilly 1991, Best 2000) than their low-redshift counterparts 
(eg. Prestage & Peacock 1988; Lilly & Prestage 1987), Best 
ct al. proposed that 3C radio galaxies at z ~ 1 and z ~ fol- 
lowed different evolutionary histories, and are seen as pow- 
erful radio sources at points within those histories at which 
they possessed similar stellar masses (~ 10 11 Mq). Under 
this interpretation, higher redshift 3C galaxies would form 
first, via merger events within their rich cluster environ- 
ments, reaching the required stellar mass by z ~ 1, before 
evolving passively to form the radio-dormant brightest clus- 
ter galaxies (BCGs) seen at z ~ 0. In contrast, the 3C radio 
galaxies at low redshift would reach the required stellar mass 
much later, via merging within their poorer cluster environ- 
ments. The argument in favour of this interpretation of the 



K — z relation was further strengthened by Best et al.'s re- 
sults suggesting that the scalelengths of z ~ 1 3C galaxies 
were significantly larger than seen in low-redshift 3C galax- 
ies, and much more compatible with the characteristic sizes 
of low-redshift BCGs. 

The evidence suggesting that 3C radio galaxies at z ~ 1 
have significantly larger scalelengths than typical of low- 
redshift 3C galaxies was re-examined using a full two- 
dimensional analysis of the Best et al. HST imaging data 
by McLure & Dunlop (2000). Comparison of the results of 
this analysis with previous results based on HST imaging 
of 3C-class radio galaxies at z ~ 0.2 (Dunlop et al. 2003; 
McLure et al. 1999) indicated that 3C galaxies at both red- 
shifts had very similar characteristic scalelengths (11 ± 2 
kpc; Ho — 50, Qo = 1). Here we can extend this comparison 
by including both the new results for the ZP5 3CRR sub- 
sample, and those for the twenty 3C galaxies included in the 
Bettoni et al. (2001) sample. 



6.1 Scalelengths 

The scalelength results for the z ~ sub-sample (Bettoni et 
al. 2001), and the z ~ 0.2 and z ~ 0.8 sub-samples (McLure 
& Dunlop 2000) were all derived by fitting classical (3 = 0.25 
de Vaucouleurs r 1 / 4 surface-brightness distributions. Con- 
sequently, in order to perform the scalelength comparison 
we also adopt the scalelength results for the ZP5 3CRR 
sub-sample obtained by restricting the fitting procedure to 
13 = 0.25. 

Under this restriction the mean scalelength of the ZP5 
3CRR sub-sample is 16.3 ± 1.6 kpc. When converted to our 
adopted cosmology the equivalent values for the z ~ 0, z ~ 
0.2 and z ~ 0.8 samples are 14.2 ± 2.2 kpc, 9.8 ± 1.6 kpc 
and 11.4 ± 2.2 kpc respectively. In conclusion, these results 
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show no indication of any trend for scalelength to vary with 
redshift in the interval 0.0 < z < 0.8. 

6.2 Evolution of the Kormendy relation 

In panel A of Fig 6 we show the Kormendy relation for all 
four 3C radio-galaxy samples. The best-fitting relation to 
the combined sample is: 

< n > e = 3.78(±0.30) logr e + 16.22(±0.33) (5) 

where all four samples have been corrected for (1 + z) 4 
surface-brightness dimming, and have had both K- 
corrections and evolutionary corrections applied as de- 
scribed in Section 5.2. 

Panel A of Fig 6 demonstrates that within the redshift 
interval 0.0 < z < 0.8 there is no obvious evidence for any 
dynamical evolution, as would be indicated by scalelength or 
surface-brightness evolution, other than that expected from 
purely passive evolution of the stellar populations. However, 
although we find no evidence for dynamical evolution of the 
class of AGN host galaxies which harbour powerful 3C-class 
radio sources in the redshift interval 0.0 < z < 0.8, dy- 
namical evolution of the individual host galaxies themselves 
within this redshift range cannot be ruled out. Panel B of 
Fig 6 shows that the redshift evolution of the mean abso- 
lute magnitudes of the four sub-samples is also consistent 
with purely passive evolution. However, it should also be 
pointed out that these results do not require that all 3C radio 
galaxies are formed in a monolithic collapse at high redshift, 
simply that the proto-galactic components which eventually 
form the host galaxies are formed reasonably co-evally, and 
that the vast majority of merger activity be completed at 
z > 1. Furthermore, neither are the results presented in this 
section inconsistent with what would be expected under the 
K — z interpretation proposed by Best et al. (1998). 

A satisfactory discrimination between the two interpre- 
tations of the K — z relation will rely on a consistent evalua- 
tion of the relative richness of the respective cluster environ- 
ments at low and high redshift. Unfortunately, much of the 
work which has previously been done on the cluster envi- 
ronments of radio galaxies, particularly at low redshift, has 
relied on statistical determinations from single-filter imag- 
ing which are known to be subject to large uncertainties 
(eg. Barr et al. 2003) . However, upon completion of our INT 
Wide Field Camera follow-up observations of the ZP5 sam- 
ple, together with our existing observations of radio-galaxy 
samples at z ~ 0.2 and z ~ 0.8, we will be able to consis- 
tently investigate the evolution in the cluster environments 
of powerful radio galaxies, by using three-colour imaging to 
successfully separate true cluster members from foreground 
and background contaminants (Jarvis et al.; McLure et al., 
in prep). 



7 COMPARISON WITH BRIGHTEST 
CLUSTER GALAXIES 

The results regarding the scalelengths and integrated lumi- 
nosities of the ZP5 sample presented in Section 5 suggest 
that the host galaxies of the ZP5 sample are directly com- 
parable with brightest cluster galaxies (BCGs). This is in 
good agreement with the results of Dunlop et al. (2003) and 
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Figure 7. Scalelength histogram for the 119 Abell brightest clus- 
ter galaxies from Graham et al. (1996). Overlaid in black is the 
scalelength distribution of the full ZP5 sample. 

McLure & Dunlop (2001b) who concluded that the hosts of 
powerful radio-loud AGN at z ~ 0.2 were directly compa- 
rable to the central BCGs of Abell class clusters. In this 
section we explore this issue further via comparison with 
the findings of the HST imaging study of z ~ 0.5 BCGs by 
Nelson et al. (2002), and the ground-based R— band surface 
photometry of 119 local Abell cluster BCGs by Graham et 
al. (1996). 

Nelson et al. (2002) investigated the characteristic sizes 
of a heterogeneous sample of 16 BCGs, drawn from both 
X-ray and optically selected cluster samples, using two- 
dimensional modelling of HST NICMOS and WFPC2 imag- 
ing with the publically available GIM2D package (Simard 
1998). The mean redshift of the Nelson et al. sample is 
< 2 >= 0.54 ±0.03, making it ideal for comparison with the 
results for the ZP5 sample. After conversion to our chosen 
cosmology, and correction for the effect of colour gradients 
(Nelson et al. 2002, Pahre 1999), we calculate the mean ef- 
fective scalelength of the Nelson et al. BCGs to be < r e >= 
14.3 ± 2.0 kpc, in excellent agreement with the scalelength 
distribution of the ZP5 sample; < r e >= 13.9 ± 0.9 kpc. 
Given that the results for both the ZP5 and Nelson et al. 
samples are derived from similar two-dimensional modelling 
of HST imaging, of samples with identical redshift distribu- 
tions, this can be taken as good evidence that the hosts of 
the ZP5 sample are indistinguishable from those of BCGs 
at z ~ 0.5. 

However, if the ZP5 scalelength results are compared 
to those determined for nearby Abell BCGs, the conclusions 
are more ambiguous. In Fig 7 we show a histogram of the 
effective scalelengths determined from one-dimensional r 1 / 4 
fits to the R— band surface-brightness profiles of 119 z ~ 
Abell cluster BCGs by Graham et al. (1996), together with 
the equivalent histogram for the ZP5 sample. The mean and 
median scalelengths determined by Graham et al. for the 
Abell BCGs are 28.7 ± 2.7 kpc and 18.5 kpc respectively 
(after conversion to our adopted cosmology). The equiva- 
lent figures for the ZP5 sample are 13.9 ± 0.9 kpc and 12.8 
kpc respectively. An application of the Kolmogorov-Smirnov 
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Figure 8. The combined Kormendy relation for the 119 Abell 
brightest cluster galaxies of Graham et al. (1996) and the full ZP5 
sample. The solid line is the best-fitting relation to the combined 
sample and has a slope of 3.01 ± 0.13. 
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Figure 9. The black-hole mass distribution for the full ZP5 sam- 
ple. The black-hole masses have been estimated from the fit- 



ted bulge luminosities via the M b ^ 
McLure & Dunlop (2002). 



M, 



bulge relation derived by 



(KS) test finds the ZP5 and Abell cluster scalelength distri- 
butions to be drawn from different parent populations at 
high significance (p = 7 x 1CT 4 ) 

In Fig 8 the Kormendy relation for the ZP5 sample 
and the Abell BCGs is shown. It is clear from Fig 8 that the 
host galaxies of the ZP5 sample are fully consistent with the 
Kormendy relation for z ~ Abell BCGs, although the ZP5 
sample, and z ~ 0.5 BCGs in general (Nelson et al. 2002), 
appear to be devoid of galaxies with scalelengths > 40 kpc. 
However, the radio galaxies comprising the ZP5 sample are 
expected to inhabit a wide range of environments (Hill & 
Lilly 1991), with previous results suggesting that the mean 
environmental richness should be comparable to low-redshift 
Abell class clusters. Indeed, if the scalelength comparison 
with the Graham et al. Abell BCGs is restricted to only 
those objects in Abell class clusters (67 of 119 objects) 
then the difference between the two distributions is only 
significant at the 2<r level (p = 0.034), in agreement with 
Hill & Lilly (1991) and McLure & Dunlop (2001b). 

In conclusion, the luminosities and scalelengths of the 
ZP5 sample show that the radio galaxies are consistent with 
the typical properties of BCGs at z ~ 0.5, which are in turn 
comparable with the central galaxies of local Abell class 
clusters. We note here that, as previously highlighted by 
Yates, Miller & Peacock (1989), 3C295 is clearly an excep- 
tional object in terms of it's cluster environment, which is 
known to be directly comparable to the richest clusters at 
low redshift (Hill & Lilly 1991). Indeed, 3C295 is a factor 
of two more optically luminous than any other member of 
the 3CRR sub-sample and, as will be seen in Section 9, also 
inhabits an extreme region of the L\^imhz — M bh plane. 



8 THE BLACK-HOLE MASS DISTRIBUTION 

As discussed in Section 5, the morphological analysis of 
the HST imaging data clearly indicates that the ZP5 radio 
galaxies are close to being classical ellipticals. Furthermore, 
it is possible to directly link the hosts of the ZP5 sample 
with z ~ ellipticals via the application of passive evo- 
lutionary corrections. Consequently, under the assumption 
that, except for passive evolution, the M b h — M bu i ge relation 
does not evolve between z = and z — 0.5 we can proceed 
to estimate the central black-hole masses of the ZP5 sam- 
ple using the correlation between black-hole mass and bulge 
luminosity observed at low redshift. 

For this calculation we adopt the M bh — M bu i ge relation 
as derived by McLure & Dunlop (2002) for a combined sam- 
ple of 72 AGN host galaxies and 18 low-redshift (z < 0.05) 
inactive ellipticals with black-hole mass measurements de- 
rived from either stellar or gas dynamics: 

log M bh / M Q = -0.50(±0.02)Af H - 2.74(±0.48) (6) 

where Mr is the absolute R-Cousins bulge luminosity (cor- 
rected to our adopted cosmology). It was demonstrated by 
McLure & Dunlop (2002) that the scatter associated with 
this relation for elliptical galaxies is entirely consistent with 
that of the M^—a relation, at only 0.3 dex (see also Erwin et 
al. 2003). Recently, Marconi & Hunt (2003) have performed 
detailed surface-brightness decompositions of 2-Micron All 
Sky Survey (2MASS) imaging of 37 low-redshift (distance < 
150 Mpc) galaxies with published black-hole mass determi- 
nations, investigating the level of scatter in the M b h — Mbuige 
relation in the near-infrared. Marconi & Hunt confirm the 
McLure & Dunlop (2002) results, extending them to demon- 
strate that the scatter associated with the M bbb — M bu i ge 
relation is ~ 0.3 dex, regardless of morphological type, pro- 
vided that accurate disc/bulge decompositions of the galaxy 
surface-brightness distributions are adopted. The distribu- 
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tion of black-hole masses for the full ZP5 sample as esti- 
mated from their bulge luminosities are shown in Fig 9. 

The geometric mean black-hole mass for the ZP5 sam- 
ple is 10 8 ' 87±0 ' 04 M Q . This is in good agreement with the 
results of Bettoni et al. (2003) who determined a geometric 
mean black-hole mass of io s - 81±0 06 M Q and io 8 s)1±0 06 M 
for 45 low-redshift radio galaxies (z £ 0.1) using the Mbh — cr 
and M bh — M bu i ge relations respectively. When combined 
with the results of previous studies at lower redshift (eg. 
Laor 2000; McLure & Dunlop 2002; McLure & Jarvis 2002; 
Dunlop et al. 2003; Bettoni et al. 2003; Barth et al. 2003) 
it is clear that truly powerful radio-loud AGN, those capa- 
ble of producing kpc-scale jets, are powered by black holes 
with masses greater than 10 Mq. In addition, given that 
radio-loud AGN are drawn from the high-mass end of the 
black-hole mass function, the results of these studies suggest 
the existence of a maximum black-hole mass of ~ 3 x 10 9 Mq 
(c.f. 3C295). This black- hole mass limit is in agreement with 
the maximum mass measured via gas or stellar dynamics in 
the local Universe (Ford et al. 1994; Tadhunter et al. 2003) 
and the limit implied from broad emission-line based black- 
hole mass estimates for > 12000 SDSS quasars at z < 2 
(McLure & Dunlop 2004). 

8.1 On the correlation between f3 and black- hole 
mass 

We note here that Graham et al. (2001) discovered a corre- 
lation between central black-hole mass and the Sersic mor- 
phological parameter f3 for a sample of 20 local galaxies with 
black-hole masses measured via gas or stellar dynamics. This 
correlation displays a comparable level of scatter (~ 0.3 dex) 
to that of the Mbh — o relation (eg. Tremaine et al. 2002) . 
Due to the fact that black-hole mass in early-type galaxies 
also correlates with bulge luminosity with comparably low 
scatter (McLure & Dunlop 2002; Marconi & Hunt 2003), 
it is expected that a correlation will exist between the (3 
parameter and host-galaxy luminosity. 

No such correlation is observed in the full ZP5 sample 
(r s = 0.20, p = 0.21). It is possible that the properties of 
two of the ZP5 sample, TOOT 0013+3459 and 3C295 (see 
Table 3), may have resulted in the determination of anoma- 
lously large (3 parameters for their luminosities. However, 
following the exclusion of these two objects the correlation 
between j3 and host luminosity is still only marginally sig- 
nificant (r s = 0.40, p = 0.012, ~ 2.4a). The principal reason 
for this discrepancy may simply be that we have a much 
smaller dynamic range in /3 compared to the Graham et al. 
(2001) study. However, we also note that the Graham et al. 
correlation between /3 and black-hole mass would imply a 
mean black- hole mass for the ZP5 sample of 1O 8 ' 4 M0, for 
our determination of < (5 >= 0.23 ±0.01. This is signifi- 
cantly different to the figure of io 8 S7±0 04 M Q determined 
from the Mbh — Mbuige relation which, using the Graham et 
al. (3 — M bh correlation, would imply < (5 >= 0.16. 



9 THE BLACK HOLE MASS - RADIO 
LUMINOSITY CONNECTION 

The question of whether or not black-hole mass is correlated 
with radio emission in AGN has been the focus of a great 



deal of attention in the literature recently (eg. Dunlop et al. 
2003; Lacy et al. 2001; Ho 2002; Woo & Urry 2002). Using 
broad emission-line based virial black-hole mass estimates 
Lacy et al. (2001) investigated the radio luminosity - black- 
hole mass relationship for a combined sample composed of 
objects from the First Bright Quasar (FBQ) and PG quasar 
surveys. The results of this study showed a strong correla- 
tion between black-hole mass and 5-GHz radio luminosity of 
the form L 5 ghz oc Af 1 ^ 4 2 over a dynamic range of some 
6 decades in radio luminosity and 4 decades in black-hole 
mass. By combining their black-hole mass estimates, derived 
from HST imaging of a sample of 33 AGN at z ~ 0.2, with 
those of the Lacy et al. study, Dunlop & McLure (2003) 
demonstrated that the 5-GHz radio luminosity of both AGN 
and quiescent nearby galaxies appears to be confined be- 
tween two radio-power limits. The lower radio-power limit 
defines the minimum radio luminosity that can be emitted 
by a given black-hole mass, and is well described by a rela- 
tion of the form L^ghz oc M^ 5 , in good agreement with the 
correlation originally observed for nearby quiescent galax- 
ies by Franceschini, Vercellone & Fabian (1998). The upper 
radio-power limit proposed by Dunlop & McLure (2003) ap- 
pears to be well described by a relation of the same func- 
tional form, offset from the lower limit by some 5 decades in 
radio luminosity. In contrast, the recent studies of Ho (2002) 
and Woo & Urry (2002) find no convincing evidence for a 
correlation between black-hole mass and radio luminosity, 
in samples comprising a range in nuclear activity from local 
quiescent galaxies up to and including powerful quasars. 

In this section we investigate whether there is a corre- 
lation between black-hole mass and radio luminosity within 
the ZP5 radio-galaxy sample. In contrast to previous stud- 
ies, the majority of which have been based on high frequency 
(5-GHz) radio luminosity, the ZP5 radio-galaxy sample al- 
lows us to test for a correlation between black-hole mass and 
extended low-frequency (151-MHz) radio luminosity. This 
distinction is potentially important given that Libimhz is 
less affected by beaming than L^ghz (eg. Jarvis & McLure 
2002) and has a close relationship to the time-averaged ki- 
netic energy of the jets (e.g. Rawlings & Saunders 1991). 

The results presented thus far are consistent with a pic- 
ture in which extended low-frequency radio luminosity scales 
roughly with host-galaxy luminosity/mass. In combination 
with the latest determination of the K — z relation (Willott 
et al. 2003) it appears that the most powerful 3C-class ra- 
dio galaxies reside in galaxies with i?— band luminosities of 
~ 4L* , while the lower- luminosity 6C and 7C-class radio 
galaxies typically inhabit hosts with luminosities of ~ 3L* 
and ~ 2L* respectively. 

At present it is unclear how the TOOT sub-sample fits 
within this picture. It can be seen from the results pre- 
sented in Table 4 that the TOOT galaxies do not appear 
to follow the rough scaling between extended radio lumi- 
nosity and host-galaxy luminosity apparent in the 3CRR, 
6CE and 7CRS sub-samples. Indeed, the mean luminosity 
of the TOOT sub-sample (3.28 ± 0.54L*) is greater than 
that of both the 7CRS and 6CE sub-samples, and is consis- 
tent with the low-redshift results of Owen & Laing (1989) 
that fat double/jet/FRI sources reside in hosts which are on 
average ~ 0.5 magnitudes brighter than those of classical 
double/FRII sources of comparable radio luminosity. How- 
ever, as was mentioned in Section 2, the TOOT sub-sample 
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Figure 10. Panel A shows 151-MHz radio luminosity plotted against black-hole mass for the full ZP5 sample. Panel B shows the 
£151 — plane for the ZP5 objects with HEG spectra. Panel C shows the L151 — plane for the ZP5 objects with classical double 
(CD) radio structures. The black-hole masses have been estimated from the bulge luminosity fits using the M bh — M bu i ge correlation 
of McLure & Dunlop (2002). The two objects highlighted with crosses are 7CRS1731+6638 and 3C295 which possess the two extreme 
values of black-hole mass and have a tendency to lead the eye. However, these two objects are included in the correlation analysis unless 
stated otherwise. 



was drawn from a preliminary version of the survey and it 
is unclear at the time of writing to what extent the current 
TOOT sub-sample is biased by the exclusion of the optically 
faintest sources. 

With this in mind, panel A of Fig 10 shows 151-MHz 
radio luminosity versus estimated black-hole mass for the 
full ZP5 sample, where the black-hole mass estimates have 
been derived via the McLure & Dunlop (2002) M bh — M bu i ge 
relation as described in Section 8. Taken as a whole, there is 
only a weak (r s = 0.35, p = 0.027, 2.2a) correlation between 
black-hole mass and radio luminosity within the ZP5 sample. 
However, as suggested previously, it can be seen from panel 
A of Fig 10 that the apparent weakness of the L151 — M b h 
correlation displayed by the full ZP5 sample is due, at least 
in part, to the inclusion of the eleven TOOT objects. 

To investigate the possible influence of radio structure 
and nuclear spectral type we have also plotted in Fig 10 
the I/151 — Mbh relation for two sub-sets of the ZP5 sam- 
ple. In panel B of Fig 10 we show the L151 — M b h rela- 
tion for those objects which display high-excitation nuclear 
spectra (HEG) only. As is suggested by the figure, this sub- 
sample displays a fairly strong correlation which is signifi- 
cant at the 99.9% level (r s = 0.62, p= 0.001, 3.0a), although 
this drops in significance to only 2.4a following the exclu- 
sion of 3C295 and 7CRS 1731+6638 (objects highlighted 
with crosses in Fig 10). In panel C of Fig 10 we show the 
L151 — M b h relation for the ZP5 objects with classical double 
(CD) radio structures only, effectively excluding the TOOT 
sub-sample. This correlation is significant at the 2.7a-level 
(r s = 0.56, p = 0.005), a result which is confirmed by the 
application of the generalized Kendall's tau test (r = 0.83, 
2.8a) , although again we note that the significance drops to 
only ~ 2a (r s = 0.43, p = 0.049) following the exclusion of 
3C295 and 7CRS 1731+6638. 

It is clear from these results that a ~ 3a correlation does 
exist between black-hole mass and low-frequency 151-MHz 
radio luminosity within the 3CRR, 6CE and 7CRS samples, 
although the scatter in radio luminosity at a given black-hole 



mass is large (~ 1.5 dex). This result is in good agreement 
with the Lacy et al. (2001) study of the L^GHz — M b h relation 
in broad-line quasars, a study which covered a much larger 
dynamic range in black-hole mass. Moreover, this result is 
also consistent with the hypothesis proposed by Dunlop & 
McLure (2003) that the radio luminosity which can be pro- 
duced by a given black-hole mass may be bounded by upper 
and lower limits separated by several orders of magnitude. 
Indeed, given that theoretical models indicate that the radio 
luminosity of AGN, in addition to black-hole mass, should 
also be related to both accretion rate and black-hole spin (eg. 
Meier 2003), it is perhaps unsurprising that the L151 — M b h 
correlation is not stronger. 

Furthermore, at present it is unclear whether the ob- 
jects which do not follow the L151 — M bh correlation, pre- 
dominantly the lower-luminosity TOOT sub-sample, do so 
because of their LEG nature, their fat double/jet radio 
structures, or alternatively, as a result of a selection effect 
which biases the TOOT sub-sample towards the inclusion 
of objects with larger host luminosities/black-hole masses. 
A more detailed analysis of the connection between the ra- 
dio properties and black-hole masses of the ZP5 sample will 
be pursued in a future paper (Mitchell et al., in prep). 



10 CONCLUSIONS 

The results of a detailed analysis of deep, high-resolution 
HST images of a sample of 41 radio galaxies which spans 
three decades in radio luminosity at z ~ 0.5 have been 
presented. The host-galaxy properties derived from two- 
dimensional image modelling have been discussed and com- 
pared to existing results in the literature. The principle con- 
clusions of this study can be summarized as follows: 

(i) The morphologies of the host galaxies of all four radio- 
galaxy sub-samples are found to be largely consistent with 
a canonical elliptical-galaxy r 1 ^ 4 surface-brightness distri- 
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bution. No trend is detected for host-galaxy morphology to 
vary significantly with radio luminosity. 

(ii) The luminosities of the host galaxies are found to be 
consistent with those of galaxies drawn from the bright end 
of the local cluster-galaxy luminosity function, in the range 
0.7 L* <L< 10L*, with a mean value of 3.2 ± 0.3L*. 

(iii) The radio-galaxy hosts are found to lie on the same 
Kormendy (/x e — r e ) relation as powerful z ~ radio galaxies, 
after allowing for the effects of passive evolution of their 
stellar populations under the assumption of a high formation 
redshift (zf or > 3). 

(iv) By combining the new results on z ~ 0.5 3CRR 
galaxies with previous results in the literature, it is found 
that there is no evidence for dynamical evolution among 
the hosts of 3C-class radio galaxies in terms of scalelength, 
luminosity and Kormendy relation in the redshift interval 
0.0 < z < 0.8. However, although no evidence is found for 
dynamical evolution among the hosts of the most powerful 
radio sources as a class, dynamical evolution of the individ- 
ual host galaxies themselves cannot be ruled out. 

(v) The scalelengths of the ZP5 host galaxies are found to 
be consistent with those of brightest cluster galaxies (BCGs) 
at the same cosmic epoch (z ~ 0.5). When compared to local 
Abell cluster BCGs the ZP5 host galaxies are found to be 
a factor of ~ 1.5 smaller on average. However, when the 
comparison is restricted to local clusters of Abell class 
only, the ZP5 and BCG scalelength distributions are found 
to be statistically indistinguishable. 

(vi) Converting the derived host-galaxy luminosities into 
black-hole mass estimates using the local M bh — M bu i ge re- 
lation predicts that the host galaxies harbour central black- 
holes in the mass range 10 8 1 M < M bh < 10 9 ' 5 M Q , with 
a geometric mean of io 8 ' 87±0 04 M . 

(vii) Significant correlations (~ 3a) are found between 
black-hole mass and extended low-frequency radio luminos- 
ity for sub-samples of the ZP5 objects which have either 
classical double (CD) radio structures or high-excitation nu- 
clear spectra (HEG). However, at present we are unable to 
determine whether the objects which do not follow the same 
correlation, predominantly the lower-luminosity TOOT sub- 
sample, do so because of their LEG spectra, fat doublc/jet 
radio structures, or alternatively, as a result of selection ef- 
fects biasing the TOOT sub-sample towards the inclusion of 
objects with higher host luminosities/black-hole masses. 
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